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I present a method for accurately calculating the pole mass of the lightest Higgs scalar boson 
in supersymmetric extensions of the Standard Model, using a mass-independent renormalization 
scheme. The Higgs scalar self-energies are approximated by supplementing the exact one-loop results 
with the second derivatives of the complete two-loop effective potential in Landau gauge. I discuss 
the dependence of this approximation on the choice of renormalization scale, and note the existence 
of particularly poor choices which fortunately can be easily identified and avoided. For typical input 
parameters, the variation in the calculated Higgs mass over a wide range of renormalization scales 
is found to be of order a few hundred MeV or less, and is significantly improved over previous 
approximations. 

PACS numbers: ll.30.Qc, ll.lO.Gh, ll.30.Pb 


Low-energy supersymmetry P| predicts the existence 
of a light Higgs scalar boson, hP , which should be acces¬ 
sible to discovery and study at the Fermilab Tevatron 
and CERN Large Hadron Collider experiments. The 
mass of hP is notoriously sensitive to radiative correc¬ 
tions. In fact, the tree-level prediction is that should 
be lighter than the boson in the minimal supersym¬ 
metric Standard Model (MSSM). It is well known that 
including one-loop corrections shows that can be heav¬ 
ier than present experimental bounds, but still leaves 
a large theoretical uncertainty, even assuming perfect 
knowledge of all input parameters. Ultimately, this sensi¬ 
tivity should become a blessing rather than a curse, since 
it means that the mass of can be a precision observ¬ 
able useful for testing particular supersymmetric models. 
This has motivated many previous efforts (for example, 
references therein) to calculate the 
higher-order corrections in various forms. 

In this letter, I will describe the calculation of the pole 
mass of hP using a method which is exact at the one- 
loop level, and includes all two-loop effects within the 
effective potential wproximation. A similar strategy has 
been employed in l^l , but neglecting two-loop effects 
involving electroweak couplings and lepton and slepton 
interactions. The complete two-loop effective potential 
has recently been given in HII^. There, I showed that 
including the previously neglected effects greatly reduces 
the renormalization-scale dependence of the minimiza¬ 
tion conditions for the Higgs vacuum expectation values 
Vu and Vd- Here I will show that there is a similar benefi¬ 
cial effect on the calculation of the mass of . Through¬ 
out, I will use the notations and conventions of M- 

The pole mass of can be calculated as follows. For a 
given choice of Lagrangian parameters specified at some 
strategically chosen renormalization scale Q in the super- 
symmetric and mass-independent DR scheme iniEi, 
one computes the effective potential in Landau gauge in 


a loop expansion: 

= + + + ( 1 ) 

One then requires that V^s is minimized [ifij to obtain Vu 
and Vd- These are Q-dependent quantities, just like the 
Lagrangian parameters, and they satisfy renormalization 
group (RG) equations which were found to two-loop or¬ 
der in M- The propagators and interactions of all of 
the fields are then obtained by diagonalizing the squared 
mass matrices, with the Higgs fields shifted by Vu and 
Vd, in the tree-level Lagrangian. This procedure ensures 
that the sum of all tadpole diagrams (including tree-level 
ones) vanishes through two-loop order. 

Then, by summing one-particle-irreducible two-point 
Feynman diagrams, one obtains the neutral Higgs self¬ 
energy matrix H^o^o (p^) at momentum p. This is a 2 x 2 
matrix (with = h^,H^) if CP violation is negligible, 
and a 4 X 4 matrix (with cjPl = /i°, G°, H°) if there is 

CP violation. For simplicity, I will assume no CP viola¬ 
tion in the following. The gauge-invariant complex pole 
mass s/jO is then defined to be the smaller eigenvalue of 

(m^Q PWhOhoijP) UhOHoijP) \ 

V n^o^o(p2) m^„+n^o^o(p2); ’ 

withp^ = syo = M’^o — iMyoVyo. The quantities m\o and 
m^o are the tree-level squared masses (without tadpole 
contributions included) as defined in section 2 of |ir)l| . 
Once the self-energy functions are known, syo can be 
found by iteration. In the following, I will quote Myo. 

In practice, the self-energies are calculated in a loop 
expansion 

n(p’) = ^n<«(p>) + ^n<»)(p») + .... (3) 

The one-loop self-energy functions n*^^)(p^) are easily 
found, but so far a complete expression for I{P\p^) is 
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lacking. However, given the n-loop contribution to the effective potential, the n-loop self-energies at = 0 are: 
_I 


V<o.o(0) 


1 fcc, -sA ( d^V‘^'^^ldvudvd\ 

2 lySa C„ ) d‘^V^^^/dvj ) 



( 4 ) 


Now, for small p^, one may reasonably approximate 
n(”)(p2) R:; n("^(0). In principle, the resulting approxi¬ 
mated pole mass suffers from two related diseases; it is 
not gauge-invariant, and as we will see it has singulari¬ 
ties (or instabilities) if evaluated at (or near) a scale Q 
at which a tree-level scalar squared mass in a loop hap¬ 
pens to vanish. However, when calculating the pole mass, 
these errors are controlled by the smallness of com¬ 

pared to the squared masses of the superpartners and 
heavy Higgs scalar bosons in loops. 

The one-loop self-energies in Landau gauge can be 
written in terms of functions: 

(5) 

SFF{m\,m\) = {m\ +m}^ - p^)Bo{m\,m}^) 

-AoimjJ - Ao(mjJ, (6) 

= 2Boim\,m\), (7) 

SsirA) = Ao{m^^), (8) 


Sv{ml) = 3Ao{ml), (9) 

Ssv{m1,ml) = (2ml - ml + 2p'^)Bo(ml,ml) 

+ (ml-p^f [Bo(ml,0) - Bo(ml,ml)] /ml 
+Ao(ml) + (ml -ml- p‘^)Ao(ml)/ml, (10) 
Svv(ml,ml) = Ao(ml)/2ml - 2Bo(ml,ml) 

+ [2(ml-p'^fBo(ml,0) - (p'^fBo(0,0) 
-(2ml - P^fBo(ml, ml)] /Mml)"^, (11) 

where 

Ho(m^) = m^[ln(m^) — 1] 

Bo(ml,ml) = 

— dx \n[xm\ -I- (1 — x)m\ — x(l — x)p'^] 

Jo 

with ln(X) = h\(X/Q^ — ie) for real X, and defined for 
complex X by Taylor expansion. Then one has: 


( 12 ) 
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fe,Z=l 


k,l^l 


+?,yl^e[A^oAl^o]SFF(t, t) + Be[k^^oAu4,o^mlS-pp(t, t)| -H 3p^|Re[A:rf,^oA:*^o]S'FF(&, b) 

+Be[kd^okd^o^mlS-pp(h, 6)| -H p2|Re[fc^^ofc*^o]S'FF(T, r) -k Re[fcd,^ofc^^o]m^S'pp(T, r)| 
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o]SrF(-^k, -^/) + ^^l^]VkNi(^9^NkNi(f>^'l''^Nk'^ -^0 | 


kj^l 


NkNi(p9-^NkN[(pO]"^jVk"^jV[^FFK 
2 2 


+3y^Re[A:^0oA:*^p] ^ Ss(tk) + ^yb^^[^d(p^^d<p9] ^ Ss{bk) + y‘^Re[k^^ok^^o] ^ Ssirt) 
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^k) 


rYg 
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fc=i 
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- K^oK^o] - xW")Ss{f) + S^j [{g^ + g'^)Sv{Z) + 2g^Sv{W)] /4 


ff" + g'- 


y^Iin[fcd,f>ofcL9 “ “ ^«0?fci^0o]>5'sy(0°, ^) 


fc=i 




fe=l 


+Re[i;„fc„0o + t)dfcd<^o]Re[z;„/c„0o + Udfe^^o] [(5^ + g'‘^f‘Svv{Z, Z) + 25"‘5'vy(VF, R")] /4. 


(14) 


The name of a particle is used to denote its squared mass when appearing as an argument of a loo p fu nction. All of 
the masses, couplings, and mixing parameters appearing here are defined explicitly in section II of [l^, except: 


—Re[kd(j,o]Re[ku^ok^^o] —Re[kd^o]Re[ky_^ak*^o]}/2V2 + (u d), 

“ {9 {b^dk^^o + Vuk^^o]kj^^+k^^++ [vdk^^o + Vukj^^o]k^^+kj^^++ SjkRe[vdkdipo + Vuky^^o]) 

~^9 [kd4,+ kd,p+ ~ kurj,+ kjj^^+]Re[vdkdti,o — Vuku,ffi]]/2\/2. (15) 


The corresponding Feynman gauge formulas are given in 
[Illli, but we need the Landau gauge results to be 
consistent with the calculation of VeS and Vu,Vd- 

The calculation now proceeds by using the above 
nR)(p^) and, as an approximation to the actual two-loop 
self-energy, the functions II^^) (0). The latter are obtained 
from eq. i) by numerically differentiating the effective 
potential appearing in M using a finite difference 
method, sampling nearby points in {vu,Vd) space. (One 
could also differentiate analytically, but the result¬ 
ing expressions are very complicated and not at all sig¬ 
nificantly more accurate.) 

Numerical results as a function of the choice of Q are 
shown in Figure ^ for the sample test model defined in 
section VI of 0- This model is defined by DR input 
parameters at a scale Qo = 640 GeV: 

g' = 0.36, g = 0.65, 53 = 1.06, 

?/t = 0.90, 2/6 = 0.13, 2/r = 0.I0, (16) 

and, in GeV, 

Ml = 150, M 2 = 280, M 3 = 800, 
at = —600, ab = —150, a,- = —40 

and, in GeV^, 

= (780)^ = (740)^ = (735)^ 

= (280)^ = (200)^ 

m|3 = (700)2, ^2^ ^ ^530)2^ ^2^ ^ (725)2, 
ml^ = (270)2, ml = (195)2, 

= -(500)2, ^2^^ ^ (270)2. (17) 


With 

fi = 504.I8II2 GeV, b = (184.22026 GeV)2, (18) 

this leads to a minimum at 

VuiQo) = 172 GeV; Vd{Qo) = 17.2 GeV. (19) 

Then the parameters of the model (including Vu,Vd) are 
run to any other scale Q using the two-loop RG equations 
of (^[l^. There, the parameters fi and b are adjusted 
to ensure that Mff is minimized; as shown in |inl| this 
readjustment is very small when the full two-loop effec¬ 
tive potential is used. Then the pole mass is found as 
described above to determine M^o, which is graphed in 
Fig. I as the solid line. Ideally, this would be indepen¬ 
dent of Q, so the fact that it is not gives an indication of 
the effects of our approximations. 

A striking feature of the graph is the presence of 
instabilities near Q = 463 GeV (where the tree-level 
squared mass of passes through 0), and Q = 568 
GeV (where the Landau-gauge tree-level squared masses 
of the Nambu-Goldstone bosons G°, pass through 0) 
0 The point is that for small tree-level scalar squared 
masses m^, the effective potential scales like 

V(2) ^ m^[cfln(m^) -b c^ln^(TO^)] -b ... (20) 

where 0^2 are constants as —> 0. Thus, while I4ff is 

well-defined and continuous in that limit, derivatives of it 
are not. (The Nambu-Goldstone bosons have C2 = 0, so 
the corresponding singularities are less severe.) These 
and nearby values of Q simply represent bad choices, 
where the approximation being made for the pole mass 
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FIG. 1: The real part M^o of the pole mass of the lightest 
Higgs boson of supersymmetry for the sample test model of 
Ref. M , as a function of the choice of renormalization scale 
Q. The solid line is the result of the calculation presented 
here. The dashed line shows the result if all effects involv¬ 
ing electroweak couplings and lepton and slepton interactions 
are removed from the two-loop contribution, corresponding to 
previous approximations. 


is invalidated by large logarithms. If it were available, 
the use of = s^o), rather than the approxima¬ 

tion n(^^(0), would eliminate the instability for choices 
of renormalization scale at which the Goldstone boson 
masses happen to vanish. [This is easily checked for 
the analogous case at one-loop order, where replacing 
by leads to similar but milder numerical 

instabilities, because of (m^)2ln(m^)/4-b ... 

for —> 0.] Therefore, one should simply be careful to 

avoid such choices for the renormalization scale 0 

For larger Q, the result for M^o is nicely stable. A 
likely good range of scale choices is 600 GeV < Q < 700 
GeV. This range includes the geometric mean of the top- 
squark masses, a traditional guess for the optimal scale 
for evaluating M^o. It also includes the scale at which 
M/jO is equal to the tree-level value m/jO, and the scale at 
which the two-loop corrections to the Goldstone boson 
masses vanish. In this range, the value of M^o calculated 
by the method described here varies by less than 100 
MeV. Even for the larger range 600 GeV < Q < 900 GeV, 
the variation of M^o is about 320 MeV. For reference, 
the precise result of the calculation at Qo = 640 GeV is 
M /,0 = II5.628 GeV in this model. 

For comparison, also shown in Figure ^ as the dashed 
line is the result which should correspond to previous ap¬ 
proximations 0,0 in which electroweak, tau, stau, and 
tau sneutrino interactions {g, g',yr,ar) are neglected in 
the two-loop part 0 Because the terms implicated in 
eq. are simply not included in this approximation, 
the instabilities of the full calculation at special values of 
Q do not appear. The more important comparison occurs 


at the better choice of larger Q as in the previous para¬ 
graph. There, the dashed-line estimate is significantly 
larger, and shows a stronger scale-dependence, than the 
calculation presented here with the complete 

I have checked that similar results are obtained in a 
wide variety of MSSM models with dimensional param¬ 
eters at or below the TeV scale, including models with 
larger and smaller tan (3 and different superpartner mass 
hierarchies and mixing angles. I find that the calculated 
MfiO is quite generally stable to within a few hundred 
MeV or less over a wide range which includes the geomet¬ 
ric mean of the top squark masses and excludes any scales 
where tree-level scalar squared masses vanish. However, 
the scale-dependence of M^o should not be confused with 
the actual theoretical error, which is probably somewhat 
larger. This is because some fraction of the neglected 
contributions is going to be scale-independent. 

To improve the situation still further, one must calcu¬ 
late the full two-loop self-energies The present 

work has shown that the effects of the electroweak cou¬ 
plings in this are certainly not negligible compared to our 
eventual ability to measure M^o at colliders. The method 
outlined here will also be a useful check on a future cal¬ 
culation of in Landau gauge, since it will have 

to coincide with the ^ 0 limit. 

The viability of any given model scenario can be tested 
by conducting global fits of Mf^o and many other observ¬ 
able masses, cross-sections, and decay rates to a set of 
underlying model parameters. If supersymmetry is part 
of our future, then the determination of M^o will play 
an important role in testing the whole structure of the 
softly-broken supersymmetric Lagrangian. 

This work was supported in part by NSF grant PHY- 
0140129. 
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